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Gradient expansion approach to nonlinear superhorizon perturbations 

Yoshiharu TANAKA and Misao SASAKI 
Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

Using the gradient expansion approach, we formulate a nonlinear cosmological perturbation theory 
on super- horizon scales valid to 0(e 2 ), where e is the expansion parameter associated with a spatial 
derivative. For simplicity, we focus on the case of a single perfect fluid, but we take into account 
not only scalar but also vector and tensor modes. We derive the general solution under the uniform- 
Hubble time-slicing. In doing so, we identify the scalar, vector and tensor degrees of freedom 
contained in the solution. We then consider the coordinate transformation to the synchronous 
gauge in order to compare our result with the previous result given in the literature. In particular, 
we find that the tensor mode is invariant to 0(e 2 ) under the coordinate transformation. 
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I. INTRODUCTION 



Recently, there have been a number of studies investigating the possibility of non-Gaussian primordial curvature 
perturbations from inflation. The cosmic microwave background (CMB) temperature anisotropy was found by COBE 
to be of order 1(T 5 [J, and a recent more accurate observation by WMAP was found to be consistent with the 
predictions of the inflationary universe, namely that the universe is spatially flat and that the primordial curvature 
perturbation is almost scale-invariant and statistically Gaussian 0, Q. Thus, it seems that the use of the linear 
cosmological perturbation theory that has been developed in the last couple of decades [HH, Q has been observationally 
justified. 

Nevertheless, the constraint obtained by WMAP on the deviation from a Gaussian distribution is still very weak, [3], 
and the next generation of observations, i.e. those to be carried by the PLANCK satellite, may be able to detect 
t— I | a deviation from a Gaussian distribution and constrain the level of such deviation by an order of magnitude better 
than the WMAP observation, with only temperature anisotropies Q, and with a combination of temperature and 
polarization (only E mode) anisotropies Q. Furthermore, it has recently been suggested that by observing 21 cm 
background anisotropies, it might be possible to constrain the deviation from a Gaussian distribution of the primodial 
curvature perturbations to an accuracy which is better than the WMAP observation by three orders of magnitude 0] . 
\ Also, theoretically, although the deviation from Gaussian behavior obtained from the standard single-field slow-roll 
inflation is too small to be detected [l^, [Til ] . many new types of inflationary models have been proposed in the 
last few years, motivated mainly by string theory and higher dimensional gravity theories, that predict rather large 
deviation s Il2t [l3| . Consequently, there have recently been many studies of the deviation from Gaussian behavior in 
inflation [li]. 

The reason that a large deviation from a Gaussian distribution has not yet been ruled out by observations lies in 
the smallness of the primordial curvature perturbation amplitude. This is demonstrated explicitly using the so-called 
Jnl parameter introduced by Komatsu and Spergel Q, 

- *i(x) + f NL ($!(*) - (*£» , (1.1) 

where 4> is the curvature perturbation on the Newton (or longitudinal) slicing and $l is the linear Gaussian part of 
the perturbation. Because <&l ~ 1CP 5 , the coefficient /jvl can be very large, even if $ is very small. 

Thus, to quantify the deviation from a Gaussian distribution and clarify its observational effect, it is important 
to develop a theory that can treat nonlinear cosmological perturbations. There are two approaches to nonlinear 
perturbations. One is the second-order perturbation theory [TEl . fl6| . which has been used by several authors to 
quantify the non-Gaussianity [lol [ill [l2l [73] ■ The other employs a gradient expansion by assuming that the spatial 
derivative is sufficiently small compared to the time derivative [l8{ . In this paper, we use the gradient expansion 
approach and formulate cosmological perturbations to full nonlinear order in its amplitude, but to second order in 
spatial gradients. 

In passing, we note that while we work with the standard metric perturbation, there exists a different approach, 
called the covariant approach [l!|. In particular, recently Langlois and Vernizzi have succeeded in formulating it in 
a form that may be useful when dealing with nonlinear perturbations not only on superhorizon scales but also on 
subhorizon scales 1201] . 
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In the (spatial) gradient expansion approach, we introduce an expansion parameter, e, and associate it with each 
spatial derivative, expand the field equations in terms of e, and solve them iteratively order by order. Originally, the 
gradient expansion was used to explore the general behavior of the spacetime near the cosmological singularity [2l|, [22j ■ 
Tomita called it the anti-Newtonian approximation [23| and used it to investigate cosmological perturbations on 



superhorizon scales 24| . The terminology, 'gradient expansion' is perhaps due to Salopek and Bond 25j. They 



adapted it to the Hamilton-Jacobi formalism and used it to investigate the nonlinear perturbations on superhorizon 
scales in slow-roll inflation. Later, Comer et al.[26| formulated the gradient expansion and discussed the form of the 
solution to 0(e ) beyond the leading order, and Deruelle and Langlois carried out a thorough study of nonlinear 
perturbations near the cosmological singularity [27| ■ 

There exists an extensive literature on the gradient expansion method. To summarize briefly, including the works 
mentioned above, gradient expansions to higher orders by iteratively solving the Einstein equations was formulated 
and studied in Ref. [HI, [H, H3, [H, [H, H3, El, [H, H3] ■ However, in most of these works, the authors either focused 
on the asymptotic behavior in the past or in the future or assumed a specific form of the solution. This obscures the 
physical degrees of freedom associated with the solutions they obtained. In particular, no clear analysis concerning 
the possible decomposition of a solution into the scalar, vector and tensor modes has been carried out. 

The Hamilton-Jacobi approach, in which one assumes the form of the generating functional to find solutions, was 
adopted in Ref. [3l|, UK, |34j. However, because of the simple form of the generating functional that the authors 
of these works assumed, the solutions they obtained were not general enough. Specifically, they ignored half of the 
physical degrees of freedom that correspond to decaying modes in the linear theory. 

Perturbations on superhorizon scales have been studied extensively with the so-called separate universe approach (35l . 
[3r| . The fact that the separate universe approach is essentially the leading-order approximation of the gradient 
expansion was demonstrated by Rigopoulos and Shellard [37| . and a more general analysis was given by Lyth, Malik 
and Sasaki (l8| . Although this leading-order approximation seems to be valid for quite a large number of inflation 
models, it can miss important points in some cases. A good such example is the case studied by Leach et al. [38| . 
They considered the linear perturbation of a single-field inflation model, which has a stage at which the slow-roll 
conditions are violated. They then showed that the 0{k 2 ) corrections to the curvature perturbation on superhorizon 
scales, where k is the comoving wavenumber of the perturbation, play a crucial role in the determination of the final 
curvature perturbation amplitude. More precisely, it is the decaying mode solution which becomes non-negligible and 
makes the 0(k 2 ) corrections important. If the 0(k 2 ) corrections are important already at linear order, one expects 
that they will significantly affect the deviation from a Gaussian behavior. In the gradient expansion, these corrections 
correspond to 0(e 2 ) beyond the leading order. Thus, it is important to include the 0(e 2 ) terms as well as the decaying 
modes in the gradient expansion. 

Here, we mention an another reason for the necessity of going beyond the leading-order approximation. While 
nothing remains unclear about the gauge dependence and/or gauge invariance in the linear perturbation theory, there 
are still some issues to be clarified in the case of nonlinear, superhorizon perturbations. Among them are the issue of 
identifying the scalar, vector and tensor modes and their gauge (or coordinate) transformation properties. Specifically, 
tensor modes which are to describe gravitational waves are not clearly seen at leading order. This is because of 
the general covariance according to which, even if there exits a gravitational wave, spacetime can be made locally 
Minkowski within a region much smaller than the wavelength of the gravitational wave. Thus, it is necessary to go to 
higher orders to clearly distinguish tensor, scalar and vector modes and to clarify the gauge dependence/invariance. 

In connection with this, we note that most authors have worked in the synchronous gauge. Of course nothing is 
particularly wrong with the synchronous gauge, but it is not the most convenient choice for the analysis of superhorizon 
scale perturbations. One apparent drawback is the existence of a gauge mode in the time coordinate. In other words, 
the synchronous time-slicing condition does not fix the slices uniquely. Thus, it is worth developing the gradient 
expansion approach in an appropriate time-slicing which fixes the time coordinate uniquely, in order to understand 
the physical meaning of solutions. 

In summary, there are several important issues in the gradient expansion approach which need more detailed 
investigations: derivation of the general solution in the sense of the Cauchy initial data; identification of the scalar, 
vector and tensor modes, if possible, in the general solution; and clarification of the gauge dependence/invariance of 
the general solution. 

In this paper, to address these issues, we formulate the gradient expansion to 0(e 2 ) beyond leading order in the 
uniform Hubble slicing. In the separate universe approach, or equivalently at leading order in the gradient expansion, 
the uniform density slicing (or equivalently uniform Hubble slicing), which fixes the time coordinate uniquely, is 
convenient [HI, [3^ |. and a nonlinear version of a gauge-invariant combination that corresponds to the curvature 
perturbation on the uniform density slicing has been constructed [l8j . This quantity is independent of the choice of 
time-slicing and conserved for adiabatic perturbations. Thus it represents a nonlinear scalar mode perturbation. This 
fact motivates us to adopt the uniform Hubble slicing. 

We employ the (3+ l)-decomposition (Hamiltonian formalism) of the Einstein equations. For simplicity, we assume 
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a perfect fluid equation of state, P = {Y — l)p, with T a constant. We then derive the general solution for all 
the variables, focusing particularly on the correspondence between the degrees of freedom appearing in the general 
solution and those in the linear theory. We find that the identification of the tensor mode in the spatial metric is 
rather arbitrary, depending on the choice of the spatial coordinates, as a reflection of the general covariance, while it 
can be unambiguously identified in the extrinsic curvature of the metric, albeit non-locally. 

Next, we consider a coordinate transformation from the uniform Hubble slicing to the synchronous slicing in order 
to compare our result with the results obtained in the synchronous gauge. The main purpose of this is to clarify the 
relation between the variables defined in each coordinate system. But this is also useful as a good way to confirm 
the computation. For definiteness, we compare our result with that obtained by Comer et al. [20|. We find that they 
are identical, provided that we choose the form of the metric in accordance with their assumption. Furthermore, we 
find the tensor mode is invariant under the change of the time-slicing. In this connection, we note that Matarrese et 
al. discussed the effect of a gauge transformation in the second-order perturbation theory in the Einstcin-dc Sitter 
universe [l5| . Comparison of their result with ours is an interesting problem left for future study. 

This paper is organized as follows. In Section II, we define the basic variables we use in our paper using the (3 + 1) 
decomposition. In Section III, we first state the basic assumptions in the gradient expansion. Then we derive the 
general solution on the uniform Hubble slicing to 0(e 2 ) beyond leading order, and identify the scalar, vector and 
tensor modes in the general solution. In Section IV, we consider the transformation from the uniform Hubble slicing 
to the synchronous slicing, and clarify the relation between the result obtained by Comer et al. (2|| and our solution. 
We conclude the paper in Section V. Some useful order-counting formulas are given in the appendix. 

II. BASIC EQUATIONS 

In the (3 + l)-decomposition, the metric is expressed as 

ds 2 — g^ydx^ dx v 

= {-a 2 + (3 k (3 k )dt 2 + 2/3 i dx i dt + j ij dx i dx j , (2.1) 

where a, (3 l ((3" 1 — J lJ /3j), and 7y are the lapse function, shift vector, and the 3-dimensional spatial metric, respectively. 
We rewrite 7y as 

lij (t,x k ) = a 2 (t)^ i (t,x k )%(t,x k ); det(7«) = l, (2.2) 

where the function a(t) is the scale factor of a fiducial homogeneous and isotropic background universe. 
The extrinsic curvature of a hypersurface of constant t is defined by 

Kij = -Vi»j = aT° tJ , (2.3) 
where n M = (—a, 0, 0, 0) is the vector unit normal to the time slices. We decompose the extrinsic curvature as 

= ^-K + i> A a 2 A l0 ■ K ee ^K l3 , (2.4) 

where A^ represents the traceless part of . The factors in front of Aij are determined so that the mixed components 
K l j take the form 

K*j = t/ 3 K + A 1 , , (2.5) 

where the indices of A^ are raised and lowered by 7 y and 7y. 
The stress-energy tensor for the fluid is 

Tfj,u = (P + P)u^u v + Pg^u, (2.6) 

where p, P and u 11 are the energy density, pressure, and the 4-velocity, respectively. We consider a perfect fluid with 
the barotropic equation of state P/p = T — 1 =const. The components of the 4-velocity are given by 

u° = [a 2 - {fi k + v k )(p k + « fc )]- 1/2 , u l = u°v\ v l ee yj'ty, 

u = -u°[a 2 -p k (p k +v k )], u i = u Q {v i +fi i ). (2.7) 



4 



We express the (3 + l)-decomposition of the energy-momentum tensor as 

E = T^n v = (p + P)(cm ) 2 - P, (2.8) 
Jj = -T^n^ = (p + P)au° Uj , (2.9) 

Sij = Tij = + P)(u°) 2 (^ + Pi^Vj + f3 3 ) + Pjij, (2.10) 

S k k = 1 kl S lk . (2.11) 
The hydrodynamic equations V M T MI/ = are written in the form [39] ] 

d t {wip 6 a 3 p 1/r ) + d k (wi/j 6 a 3 p 1/r v k ) = 0, (2.12) 
d t {w^ e a 3 (p + P) Uj } + d k {wt/j e a 3 (p + P)v k Uj } 

= -atp 6 a 3 d 3 P + wip 6 a 3 (p + P){~au°d 3 a + u k d.jf3 k - y^<9j7 fc '}, (2.13) 
where w = au°, and we have 

.fc _ nk , 

iP 4 a 2 u° 
from Eq. (|2T7]) . 

In the (3 + l)-formalism of the Einstein equations, the dynamical variables are jij and Kij. The (n,n) and (n,j) 
components of the Einstein equations give the Hamiltonian and momentum constraint equations, respectively, while 
the components gives the evolution equations for K^. The evolution equations for 7^ are given by the definitions 
of the extrinsic curvature (|2.3|) . 

In the present case, the Hamiltonian and momentum constraints are 

R - Aij A j + \k 2 = WirGE , (2.15) 

2 
3 

The evolution equations for 7^ are given by 



DiA^ - -DjK = 8nGJj . (2.16) 



(d t - P k d k )^ + j;{-aK + 8 k p k }, (2.17) 

2a 6 

(d t - P k d k )^ = -2aAa + i %k d 3 f3 k + ^ Jk d t /3 k - pijd k f3 k , (2.18) 



where the dot represents d/dt. The evolution equations for are given by 



[d t - /3 k d k )K = aiAijAv + ^K 2 ) - D k D k a + AvGa{E + S k k ), (2.19) 



(d t -(3 k d k )A l3 = J—[ a (R t3 -l^R)-( Dl D 3 a-^D k D k a)] (2.20) 



+a(KA ij - 2A lk A k 3 ) + A lk dtf k + A jk d t f3 k 
2 A f) R k 87rGa /c Hi ok x 



where R^ is the Ricci tensor of the metric 7«, with R = r ) lJ Rij, and Dj is the covariant derivative with respect to 

We define the local Hubble parameter as 1/3 of the expansion of the unit normal vector n M , which is equal to —1/3 
of the trace of the extrinsic curvature in our convention: 

3a d t ib Di/3 l 

3H=-K = — + 6— — . (2.21) 

aa aip a 

In the following, we adopt the uniform Hubble slicing. For this slicing, we have 

H(t) = ~ , (2.22) 



and Eq. (f2~2T]) implies 



Q = 1 + W^w- (2 ' 23) 
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III. GRADIENT EXPANSION 



A. Assumptions 

We consider nonlinear superhorizon perturbations in the approach called the anti- Newtonian approximation [2^ |. 
the spatial gradient expansion [25| , or the long wavelength approximation [3^ | . In this approach, we assume that the 
characteristic length scale L of inhomogeneities is always much larger than the Hubble horizon scale, i.e. L 3> H^ 1 ~ t. 
We introduce a small parameter e, and assume that L is 0(l/e). This assumption is equivalent to assuming that the 
magnitude of spatial gradients of the quantities is given by diip — ip x O(e), did = a x 0(e), and so on. In the limit 
L — > oo, i.e., e — > 0, the universe looks locally like a FLRW spacetime, where 'locally' means as seen on the scale of the 
Hubble horizon volume. It is noted that physical quantities which are approximately homogeneous on each Hubble 
horizon scale can vary nonlinear ly on very large scales. Thus this approach at leading order is called the separate 
universe hypothesis (3f| Hi| . 

The local homogeneity and isotropy imply that anisotropic quantities, i.e., (3 Z and v l should be at least O(e). Also, 
for the spatial metric to be locally homogeneous and isotropic, 7^ should be time independent in the limit e — > [l8| . 
This implies d t %j = 0(e). As discussed in § IIII1 B. allowing it to be O(e) introduces a decaying mode at this order, 
which affects all the quantities at 0(e 2 ) and makes the analysis too complicated. However, as discussed there, we have 
dtjij = 0(e N ) with TV > 3 in the case of an inflationary universe. Therefore, we expect that requiring a bit tighter 
condition, dtjij = 0(e 2 ), will not make our analysis meaningless. In other words, the assumption dtjij — 0(e 2 ) 
should be sufficiently general for our analysis to be applicable to a sufficiently wide class of cosmological situations. 
Thus our basic assumptions are 

I? = 0(e), v* = 0(e), d t %=0{e 2 ). (3.1) 

We note that from Eq. (|2.f 8[) . the assumptions p3 l — O(e) and dtjij = 0(e 2 ) imply A %j = 0(e 2 ). 
Next, we introduce 

S = p/p -l, (3.2) 
where po is a fiducial homogeneous part of p, which is defined by 

Po *) = • (3.3) 

07TG 

Here, H(t) is given by Eq. (|2.22p . Then, using the equations presented in the previous section, we evaluate the orders 
of magnitude of the basic variables as [391 ] 

iP = 0(l), x = a-l = 0(e 2 ), w = au° = I + 0(e 6 ) , 
d t % = 0(e 2 ) , dti> = 0(e 2 ) , A l3 = 0(e 2 ) , 

5 = 0(e 2 ), u 2 = 0(e 3 ), (f +v l = 0(e 3 ) . (3.4) 

Except for dtjij = 0(e 2 ), the orders of magnitude of the above quantities can all be derived from the Einstein 
equations. An explicit demonstration of the above evaluations is given in the appendix. 

In particular, the fact that ip is constant in time to O(e) follows from the perfect fluid assumption, P = P(p). This 
corresponds to the conservation of the adiabatic curvature perturbation in the linear theory. Hence, we may regard 
-0 as a nonlinear generalization of the curvature perturbation, although the actual spatial curvature R is no longer 
described solely in terms of ip in the nonlinear case. 

We also note that the curvature perturbation is defined in several different ways in the literature. The correspon- 
dence between ip used here and the conventions used in two other papers are 

■p 2 = exp^uvis] = exp[( Ma i] , (3.5) 
where "0lms and Civiai are the quantities used by Lyth et al. [l8| and Maldacena [ll] ]. respectively. 



B. The leading-order solution 



First we derive the leading-order solution for our basic variables. We note, however, that this does not mean that 
we only solve the Einstein equations at lowest order in the gradient expansion. For example, from Eq. (13. 4p . we 
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have Aij = 0(e 2 ). Hence we may ignore it at lowest order in the gradient expansion, while here we do solve for the 
leading-order part of Ay . 

Taking the 0(e°) part of Eqs. (|2T5)) and (|2~1"9|) , we have 

~K 2 = 16^Gp , (3.6) 

8 t K = ix 2 + 47rG(3r-2)p . (3.7) 

Integrating these equations, we obtain 

a(t) = a f t^, (3.8) 

^ = 6^gW' (3 - 9) 

where aj is a constant whose normalization is arbitrary. 

Substituting the order of magnitude evaluations of the variables shown in Eq. (|3.4p into the hydrodynamic equations 
given in §2, we find 

Us + ^ + DkV * = (e% (3.10) 
d t {^p a 3 T Uj } = -TP 6 p a 3 [Td jX + (T - 1)8,6} + 0(e 5 ) . (3.11) 
The Hamiltonian and momentum constraint equations give 

j ij DiDjtp = \i a Rki^ ~ 2TrGiP 5 a 2 Po 6 + 0(e 4 ), (3.12) 

o 

Di^Aij) = 8irGT poUiip e + 0(e 5 ), (3.13) 

where -Ry is the Ricci tensor with respect to 7y , is the covariant derivative with respect to 7y . 
The evolution equations for the spatial metric give 

6 -^-3 a - X = D k ^, (3.14) 

(d t - (3 k d k )^ 3 = -2Ay + j lk d 3 (3 k + jjkd^ - ^ l3 d k (3 k + 0(e 4 ) , (3.15) 



while the evolution equations for the extrinsic curvature give 



a ~ 1 7y 

V 2 X = ±TrGp Q a 2 {(3T - 2)6 + 3T X } + 0(e 4 ) , (3.17) 



dtAii + 3- Ay = [i2y - fR] + 0(e 4 ) , (3.16) 



V^^^V^). (3-18) 



where 



We now solve the above equations to derive the leading-order general solution for each variable. From Eq. 
we have 

^ = {0) L(x l ) + O(e 2 ), (3.19) 

where ( )L is an arbitrary function of the spatial coordinates and is O(e ). Here and in what follows, the prefix (n) 
means that the quantity is 0(e n ). 

From Eq. (|3~TU|) and Eq. (pTL4"|) , we obtain 

d t 6 + 3T- x = 0(e 4 ). (3.20) 
a 
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Since the left-hand side of Eq. (|3. 1T[) is 0(e 4 ), the fact that both 5 and \ are 0(e 2 ) gives 



From Eqs. (f3T20|) and ([3T21]) . we have 

Integrating this equation yields 
Then, Eq. (j3~2Tj) gives 



W — 2 

X= ^<5 + 0(e 4 ). (3.21) 



5t<5-(3r-2)-5 = 0(e 4 ). (3.22) 
a 



5 = (2) Q(x l )t 2 -^ +0( e 4 ). (3.23) 



X = -^ F ^(2)Q(^)^*+0( £ 4 ). (3.24) 



Next, we seek the general solution for m. First, using Eq. (|3.21|) , Eq. (|3.11[) becomes 



9 t {p a 3 r % } = + (e 5 ). (3.25) 



Solving this equation, we obtain 

d ^ Qix %^+ {3)Cj (xy 2 -r, (3.26) 



J 3r + 2 

where (3)6^ (x 1 ) is an arbitrary vector of the spatial coordinates. 

To find the leading-order solution for A^, we need to evaluate the spatial Ricci tensor R^. For this purpose, we 
first consider the form of the spatial metric. Noting the relation dt^ij = 0(e 2 ), we define 

H^t, x l ) = 7y - (0 )h 3 (x l ) = 0(e 2 ), (3.27) 

where (o)fij is an arbitrary function of the spatial coordinates. Namely, Hij(t, x r ) is the 0(e 2 ) part of 7^. Then the 
Ricci tensor of the spatial metric ip^Jij can be decomposed as 

Rij = ^ + R% , (3.28) 

where is the Ricci tensor with respect to 7^ , and 

2 2 6 2 

R% = --DiDjiff - -7y At/. + —D^D^ - —^D^D^ , (3.29) 

where A is the Laplacian with respect to 7^. Using Eqs. (|3.27[) and (|3.19[) . we have 

Rij = ( 2 )^+0(e 4 ), 

4 = (2)R^+0(e 4 ), (3.30) 
where (2) Rij is the Ricci tensor with respect to (o)/ij, and 

2 - - 2 

-j D i D j(a)L 

(0)^ (o)-k 



(2)Rfj = rA5j( )L - (0) /.yA( )L (3.31) 



6 - - 2 - fc 

H T2 D i(a)LDj(p)L —^f^Dk^LD (0 )L, 

(0)^ (0)^ 

where is the covariant derivative with respect to ( ) /y , and A is the Laplacian with respect to ( ) /y . 

Now, we seek the general solution for Aij. Using Eqs. (|3.28[1 and (|3.30p . the evolution equation for Aij, Eq. (|3.16D . 
becomes 

d t Aij + ^Aj = (2) y } + 0(e 4 ), (3.32) 
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where 



(2)Fij 



(0) 



L 4 



(2)Rij + (1)Rfj ~ ^7^- (0)f kl {(2)Rkl + (2)Rkl) 



Equation (|3.32p can be immediately integrated, and we obtain 



Aij — 



a(t')dt'+ i2)C % {xk) +0(e% 



(3.33) 



(3.34) 



where (2)Cy is an arbitrary symmetric, traceless tensor of the spatial coordinates. Thus we obtain 



A, 



ST 



a 2 f (3r 



(2) F ij 



, h\ 1— A. (2)C-ij(x k ) _£ 
(X )t 3r + U ^ if r 



(3.35) 



Here we note that, because the last term proportional to (2)CV/ a 3 is a homogeneous solution to Eq. (|3.32[) , it may 
be 0(e) in principle. However, allowing it to be 0(e) alters the 0(e 2 ) part of the Hamiltonian constraint (|3 . 1 2[) 
substantially, because the quadratic term AijA lJ would no l ong er be negligible. This situation is essentially the same 
as that in the case of homogeneous, anisotropic cosmologies [22l l24l l27j| . in which the square of the shear (the traceless 
part of the extrinsic curvature) acts as an effective energy density proportional to a~ 6 . To avoid this complication, 
for the sake of simplicity, we have assumed dtjij — 0(e 2 ), which implies A^ = 0(e 2 ). We have not investigated 
how restrictive this additional assumption is. However, in the case of an accelerating universe, < V < 2/3, which 
would mimic an inflationary universe, this decaying mode (oc a -3 ) is expected to be of order H(k/(aH)) N , where 
N = 3(2 - T)/(2 - 3r), so that H(k/(aH)) N oc a~ 3 . Since N > 3 for < T < 2/3, it should be higher than third 
order in e. Thus, the assumption A^ = 0(e 2 ) should be valid. 

Up to now we have not considered the Hamiltonian and momentum constraint equations. The constraint equations 
will relate the arbitrary spatial quantities (a)L, (o)fij, <2)Qi (2)Cjj and (3)Cj- 

The 0(e 2 ) part of Hamiltonian constraint (j3 . 12[) yields 



(2) 



Q 



3r 2 



t 



-(o)f 3 DiDj^L + -(o)f kl (2)Rki(o)L 



+ 0(e% 



a f(0) Ll 

-X-TI (0)f l [(2)R L u + (2)Rm] + 0(6 4 ) = R [L 4 /] + 0(e 4 ) 



f 



(3.36) 



where R [L 4 f~\ is the Ricci scalar of the metric ( )i 4 (o)/ij- Thus (2)Q is not arbitrary, but is expressed in terms of 
(o)fij and ( )L. 

The 0(e 3 ) part of the momentum constraint (|3.13p yields 



LP 



(0) 



= 8TrGTp u m L 6 + 0{e 5 



(3.37) 



Inserting Eqs. (|3.35p and (|3.26p into this equation, we find 



d i(2 )Q = 9r 2 

3r + 2 4aJ(3r + 2) (0) i 6 



D j [(o)i 6 (2)^] , 



(3)Ci 



3r 



4a 3 (0) i 6 



D 3 [( )i 6 (2)Cy] 



(3.38) 
(3.39) 



The latter equation implies that (3)Cj is not arbitrary but expressed in terms of (o)L, (o)/ij and (2)Cij, while the 
former equation is found to be consistent with Eq. (|3.36|) . This consistency is a result of the Bianchi identities. 

Before closing this subsection, we note the following fact. We have considered Uj in the above, but we have not 
considered the 3- velocity v l and the shift vector (3 l individually. As seen from Eq. (|2.14| . we have i/ + (f = vPj^Uj. 
Thus, once we fix the shift vector under the assumption (3 l = 0(e), v l is immediately determined. The reason that 
we have not bothered to do this is that it was unnecessary to fix the shift vector for the derivation of all the other 
quantities in order to derive their leading-order terms. Geometrically, Uj does not depend on the choice of the shift 
vector because it describes the components projected on to the hypersurfaces of constant t. The same is true for all 
the other variables, except for v l . This shift vector independence is probably a result of the localness of the gradient 
expansion. 
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C. The solution to 0(e 2 ) in the gradient expansion 

Now we consider the general solution valid to 0(e 2 ) in the gradient expansion. As we have seen in the previous 
subsection, among the basic variables we have introduced, the only quantities whose leading-order terms are lower 
than second order in e are i\) and 7^, that is, the spatial metric. Hence, we have to evaluate the next order terms of 
these variables. 

To find the next order solutions for ip and 7^ , we now have to specify the shift vector. For simplicity, we choose 
(3 l to be 0(e 3 ). Under this assumption, the solutions are easily derived. We note that this choice includes the case of 
comoving coordinates, v % — 0, as seen from Eq. (|3.4[) . 

First, we consider t/j. From Eq. (|3.14|) . we have 

^-3^=0(e 4 ). (3.40) 
tp a 

Substituting the solution (|3.24|) for \i this equation is readily integrated, and we obtain 

* - mi - + c^o*') + ( 3 - 41 ) 

where ( 2 )L is an arbitrary function of the spatial coordinates. Without loss of generality, we can absorb (2)L into 
(q)L. Thus we obtain the final expression for r/> as 

^ = (0) X(,V (0) ^ 2)0( "% 2 -^ + Q(e 4 ). (3.42) 
Turning to 7^, from Eq. (|3.15[) . we have 

dtjij = -2Aij +0(e 4 ). (3.43) 
Substituting the solution (|3.34[) for Ay, this equation gives 

4 dt f* , , , „ f* dt 



Ha = -2 w F tj I J a(t)dt 2 (a) C« J ^ + 0(e 4 ), (3.44) 



The integrals are easily evaluated, and we find 

_ (2)Fjj 9r 2 2 __^ 2r (2) c u - 2 k 4 

4 9r2 _ 4 * 3r - a 2 (r -2)* + ^ ) + 0(e)i (3 ' 45) 

where (2)/ij is an arbitrary function of the spatial coordinates. The final expression for 7^ is given by 



k { 2)F^ k ) 9r 2 , 2 _^ 2r (2) c^) , 

7« = (o)/y(z ) -5 a^T" 7* * ~ „2/ r _ ^ * r + °( e ) 



a 2 , 9r 2 -4" Oy(T-2) 



(3.46) 



where we have absorbed (2)/y into (o)/y- Here, it is important to note that (o)/ij is not completely arbitrary, but its 
determinant must be unity: det((o)/y) = 1 + 0(e 4 ). 



D. Summary of the general solution 

We now have the general solutions valid through 0(e 2 ) for all the physical quantities on the uniform Hubble slicing, 
under the condition /3 4 = 0(e 3 ). We list them below: 

a = 1+X = 1-(2)Q^=^* 2 -^, (3.47) 

V> = (o)L (l (2) Q ^ 2 ^) = (0)L (l [L 4 /] t 2 ^ J , (3-48) 
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h = (^ fl ^ +2) ^+(^^-*. (3-50) 

^ = ( 2) Qi 2 ^, (3.51) 

qp2 op 
«< = (o)^ 6 ^ [(0)i 6 ( 2 )^] 4fl 2 (3r + 2) ^ + (o)i- 6 ^' [(o)i 6 (2)^] t 2 -r , (3.52) 

where the freely specifiable functions are (o)L, (o)/ij and ( 2 )CV/- The functions (2)Fij an d (2)Q are given by Eqs. (|3.33|) 
and (|3.36p , respectively, as functions of r Q \L and (o)/ij- Note that ( )L 4 ( 2 )-Fij is the traceless part of Rij [L 4 f] , i.e., 
the traceless part of the Ricci tensor of the metric (o)£ (o)/tn as seen from its definition (]3 .33|) . 

To clarify the physical role of these freely specifiable functions, let us first count the degrees of freedom. Because 
the determinant of (o)/ij is unity, it has 5 degrees of freedom, and because ( 2 )CV; is traceless, it also has 5 degrees of 
freedom. Thus, together with the degree of freedom of (q)L, the total number is 1 + 5 + 5 = 11, while the number of 
the true physical degrees of freedom is 4 + 2 + 2 = 8, where 4 are for the fluid (1 from the density and 3 from the 
3- velocity) and 2 + 2 are for the gravitational waves (2 from the metric and 2 from the extrinsic curvature). This 
implies that there remain 3 gauge degrees of freedom. It is easy to understand that these 3 gauge degrees come from 
spatial covariance, that is, from the gauge freedom of purely spatial coordinate transformations, x % — > x 1 = / l (x J ). 
Thus we may regard m) fij as containing these 3 gauge degrees of freedom. 

To summarize, the degrees of freedom contained in the freely specifiable functions can be interpreted as 

(0 )L • • • 1 = 1 (density) , 
(0)/« ••• 5 = 3 (gauge) + 2 (GWs) , 

(2) Cy ••• 5 = 3 (velocity) + 2 (GWs) . (3.53) 

To fix the gauge completely, we therefore must impose 3 spatial gauge conditions on (o)/ij- If (o)/y were not a metric, 
we could impose the transversality (divergence-free) condition on it. However, this cannot be done, because (o)/ij is 
the metric and hence, any covariant derivative of it vanishes identically. (This is the very nature of any metric tensor.) 

However, in the case of linear perturbation theory, or even in the case of a higher-order perturbation theory, this 
difficulty disappears because of the presence of the background metric. In such a case, we may impose the transversality 
condition on the full (nonlinear) metric with respect to the background metric to single out the gravitational wave 
degrees of freedom. In fact, the most commonly used conditions are those realized by imposing transversality with 
respect to the flat Cartesian metric, 

<5 fc ^ fc(0 )/ y - =0>' (O )/y =0. (3.54) 

At this stage, it is useful to clarify the correspondence between the degrees of freedom counted in Eq. (|3.53p and 
those of the linear theory. From the time dependence associated with (p)L, (p)fij a nd ( 2 )CV,-, we can identify ( )L with 
the growing adiabatic mode of density perturbations, the 2 degrees in (o)/y with the growing modes of gravitational 
waves, the 3 degrees in ( 2 )Cy with 1 decaying scalar and 2 vector (vorticity) modes, and the remaining 2 in (2)Cy 
with the decaying gravitational wave modes. 

Now, returning to the full nonlinear theory, it is generally impossible to identify uniquely the gravitational wave 
degrees of freedom in the metric, because of general covariance. Nevertheless, if we focus on the extrinsic curvature 
Kij , since its transverse-traceless part may be defined unambiguously, we may be able to single out the gravitational 
wave (tensor) degrees of freedom. 

To show that this is indeed the case, let us consider the momentum constraint (|3.37[) . 



D J 



(o)L 6 Aij 



= 8irGTp a Ut {a) L 6 + (9(e 5 ). (3.55) 



If we require the tensor modes to be transverse-traceless, an apparent candidate is the part of ^L 6 Aij that does not 
contribute to the right-hand side of the momentum constraint. In other words, we identify the transverse-traceless 
(TT) part of ^L 6 Aij with respect to the metric (o)fij as the tensor modes. Then, the question is if the TT part can 
be uniquely determined. Fortunately, it is known that this is possible [40(. Here we recapitulate the result briefly. 
For any symmetric traceless second rank tensor Xij (with respect to 7^ in our case), we have the decomposition 

X tj = DiWj + b l W i — '^ji J D k Wu + Xj^ 

= (AWh+X^, (3.56) 
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where W% satisfies the second-order differential equation 

[O W)i = —&{KW)ij = -Di(Xij) . (3.57) 

York showed that the solution can be determined uniquely [4(ij |. He then further showed that this decomposition is 
conformally invariant. That is, for 7^ = <j> 4 ^fij, we have = <fi~ 2 Xij, and therefore as the TT part, 

Xj T = r 2 Xj? . (3.58) 

Thus, identifying Xy with ^L e Aij, we immediately see that the transverse parts of ( )L 6 ^-Fy and (o)-^ 6 (2)Cy 
represent the gravitational wave modes. Furthermore, from their time dependence, it is clear that the former is the 
growing mode while the latter is the decaying mode, as in the case of linear theory. 

Here, it may be worth noting that the TT part of (o)L 6 {2)Pij exists even if the metric (o)/y is trivial, say, (o)/y = 
The non-vanishing TT part in this case is due to nonlinear interactions of the scalar adiabatic growing mode, as can 
be seen from Eq. (|3.33[) . 

Finally, let us comment on the vector mode. It is contained only in (2)Cy , as can be seen by considering the vorticity 
conservation law, which is valid for any perfect fluid [4l|. It can be shown that the time dependence i 2 ~r oc a 3 ^ r_1 ' 
associated with (2)Cy is precisely that resulting from the vorticity conservation. 



IV. COORDINATE TRANSFORMATION TO SYNCHRONOUS COORDINATES 

Comer et al. [26] obtained the metric to 0(e 2 ) in the synchronous coordinates. However, the result they obtained 
[Eq. (2.14) in their paper] has a rather different form from ours, partly because of the difference in the choice of 
coordinates and partly because of the fact that their result is not the general solution. Here, we seek a coordinate 
transformation from our coordinates to the synchronous coordinates and derive an explicit relation between their 
solution and ours. 

Let us denote the coordinates in the uniform Hubble slicing by {x^} and those in the synchronous gauge by 
and the corresponding spacetime metrics by g^(x) and g^ v {x), respectively. The solution obtained by Comer et al. 
is 



9ij 



= t 3r 



3T- 



4(9r - 4) 

9r 2 



9r 2 



t 2 -^ + 7 ' 2 t- 1 +j 2 )h l:j R[h} 



(4.1) 



where 



(0) L 4 (o) fij 



(4.2) 



and Rij[h] and R[h] are the Ricci tensor and Ricci scalar, respectively, of the metric hij, and P2, P2, 72 and 79 are 
arbitrary constants. Here and in the rest of this section, we set a t = 1, in accordance with the convention of [261 ] . We 
note that the terms containing /? 2 and 72 can be absorbed into fty , and hence we may effectively set these constants 
to zero without loss of generality. Thus the metric we consider is 



3T 2 



4(9r - 4) 

9r 2 . 

n f 

9r 2 -4 



hijR[h] 



R lj [h] --^R[h] 



(4.3) 



It should be noted that the term 7 2 ^ _1 in the first line of the above is a gauge mode that exists in the synchronous 
gauge, corresponding to the shift of the time slice t — *■ t + 8t(x l ). 
Now consider the coordinate transformation 



(4.4) 
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Under this transformation, the metric is transformed as 



, . dx a dx 13 _ ._. 

For g^ v being the metric in the synchronous coordinates, that is, for <?oo = — 1 and <?o 

df dp dp dp 



(4.5) 



0, we have 



5oo (x) 
9ko(x) 
9ki{x) 



-(l + 2 X )+0(e 4 ) 



dt dt dt dt 



9ij(x) , 



fl n, ^ df df dp dp _ 



df df dp dp 
dx k dx 1 dx k dx 1 



9ij(x). 



(4.6) 



Since the two metrics should coincide in the limit e — ► 0, and the non-trivial corrections are 0(e 2 ) in both metrics, 
we may choose the coordinate transformation to be of the form 



(4.7) 



where F M = 0(e 2 ). Here we have assumed that the spatial coordinates also coincide in the limit e — > 0, although 
there always exist 3 purely spatial coordinate degrees of freedom, as discussed near Eq. (|3.54[) . Inserting the above 
into the first two relations in Eq. (|4.6[) . we readily find 



dF° 
~dt 



= X. 



and 



while the last in Eq. 



i 2 h 



dF^ 



dF° dF i _ _ 2 tf dF° 

dx 1 dt dxi 



(4.8) 



(4.9) 



gives 
9ij{x) 



gij(x) + 0(e 3 ) = g i:i {x) + F°(x)d t g lJ (x) + 0(e 3 ) 
g lJ (x) + 2H(t)F Q (x)g lJ (x)+0(e 3 ). 



(4.10) 



Here it is important to note that this coordinate transformation of the spatial metric is a conformal transformation 
to 0(e 2 ). Thus it leaves the unit determinant part of the metric, jij, and its time derivative (or the traceless part of 
the extrinsic curvature), Aij, unchanged. Recalling that the TT part of Ay with respect to the metric 7y is uniquely 
determined and is invariant under a conformal transformation of the metric, up to conformal rescaling of the TT 
part, it seems reasonable to conjecture that the tensor modes are time-slice independent to 0(e 2 ) for a wide class 
of time-slicings, though the above result is valid only for the transformation between the synchronous and uniform 
Hubble slices, in the strict sense. 

Equations (|4.8[) and (|4.9[) determine the coordinate transformation. Integrating Eq. (|4.8p . we find 



F° = / X dt = - 



3r - 2 



9r 



(2)Qt 3 



(4.11) 



3r 2 (3r - 2) 

8(9T-4) 



R[h]t 



{2)E{x i ), 



where (2)E(x l ) is a function of the spatial coordinates (determined below), and the second equality follows from 
Eq. (|3~3^)l . Substituting this into F° in Eq. P~TU]) and noting Eq. (|3~3"]) . we find the metric in the uniform Hubble 
slicing to be 



9ij(x) 



t3T 



hij (l - ~R[h] t 2 -^ + ^R[h] + ^(2)^} t- 1 



9T 2 



9T 2 



t 2 ~sr +f3' 2 p 



R %3 [h] 



-R[h] 



(4.12) 
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Apparently, the terms in the curly brackets proportional to i -1 should vanish for uniform Hubble slicing. This 
condition determines the function n\E{x % ). In passing, we note that the solution obtained by Comer et al. assumes 
the gauge mode to be proportional to R[h], but it can be an arbitrary function of the spatial coordinates in general. 

Now, from Eqs. (|3.48[) and (|3.49p , and noting that — tsrip 4 ^^ and /ly = (o)^ 4 (o)/iji our solution can be 
expressed as 

9ij {x) = mh ik (l- ^-R[h]t 2 -^) (4.13) 



where from Eq. (|3.33|) . we have 
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A fe 9r 2 k 2-4, . x _ 



pjF^fclftl-^j. (4.14) 
Comparing this with Eq. (|4.12[) . we see that these two metrics are identical, provided that we choose 

2r 



r-2 



(2)C fc i=# (a )FV (4.15) 



Thus, we have clarified the relation between our solution and the solution obtained by Comer et al. [2g|. As is clear 
from this equation, their solution corresponds to a special choice of (2)C»j> similar to the choice of their gauge mode 
mentioned above. Although the choice of (2)Cy is unimportant in the present case of a perfect fluid because it contains 
only decaying modes, such a special choice is not only required but also renders the physical meaning of (2)Czj unclear. 



V. CONCLUSION 



In this paper, focusing on the simple case of a perfect fluid with constant adiabatic index, we have considered 
nonlinear perturbations on superhorizon scales in the context of the (spatial) gradient expansion. We have adopted 
the uniform Hubble slicing and derived the general solution valid to second order in spatial gradients. This general 
solution is found to contain 11 arbitrary functions of the spatial coordinates. We have identified the physical meanings 
of all these functions. In particular, the relations between these degrees of freedom and those in the linear perturbation 
theory have been clarified. They are as follows: 

1 • • • growing scalar (curvature) perturbation , 

1 • ■ • decaying scalar perturbation , 

2 ■ • • decaying vector (vorticity) perturbations , 
2 • • • growing tensor perturbations , 

2 • • • decaying tensor perturbations , 

3 • ■ • spatial gauge degrees of freedom . (5.1) 

In particular, we have found that the tensor modes contained in the extrinsic curvature (the time derivative of the 
metric), which can be identified using York's conformally invariant decomposition method [4fl |. are determined non- 
locally. Furthermore, we have found that the tensor modes can be generated by nonlinear interactions of the scalar 
growing mode, even if there exist no tensor growing modes, i.e., (o)/jj = Obtaining a detailed understanding of 
this result and its implications is left for future studies. 

Then we compared our result with that obtained by Comer et al. in the synchronous gauge [26j by finding the 
coordinate transformation from the uniform Hubble slicing to the synchronous (proper time) time-slicing. In carrying 
out the coordinate transformation, we found the time-slice independence of the tensor modes to 0(e 2 ). Then, we 
found that the solution obtained by Comer et al. corresponds to a special choice of 5 = 1 + 2 + 2 decaying modes 
(i.e., 1 scalar, 2 vector and 2 tensor decaying modes) of the 11 arbitrary functions mentioned above. 

Although these decaying modes are unimportant in the present case of a perfect fluid, they may not be negligible in 
the case of a scalar field. In particular, as mentioned in the introduction, the decaying modes may play an important 
role, even in the case of a single-component scalar field, if the conventional slow-roll condition is violated. In this 
respect, one of the advantages of the framework presented in this paper is that it can take account of all the physical 
degrees of freedom. 
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Finally, we note that our result contains not only nonlinear scalar modes but also the nonlinear coupling of scalar 
and tensor modes. Thus, although the present paper is limited to the case of a simple perfect fluid, if the framework 
we have developed is extended to the case of a (multi-component) scalar field, the resulting formalism will be a 
powerful tool to investigate various non-linear, non-local effects in inflationary cosmology, especially the deviaton 
from a Gaussian behavior of curvature perturbations in inflation. We will attempt to develop such an extension in 
the near future. 
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APPENDIX A: ORDER COUNTING 



Here we demonstrate explicitly the order counting of the basic variables given in Eq. (13. 4p . This is essentially a 
recapitulation of the analysis given in Ref. (39j . 
In the (3 + l)-decomposition, u° is represented as 

u° = [a 2 -(f3 k + v k )(f3 k + v k )]-^ 2 . 

Then, from the first two basic assumptions in (|3.1|) , i.e., if = 0(e), (3 l — 0(e), we have 

u = - + O{e 2 ). (Al) 
a 

This equation and Eq. (|2.8j) yield 

E=(p + P)(au°) 2 -P = p + 0(e 2 ). (A2) 
Thus, from the Hamiltonian constraint (|2.15[) . we obtain 

H 2 = ^ + 0(e 2 ), 

= gZg^l±j)+Q(e 2 ). (A3) 
3 

Now, since H is uniform, by definition, on uniform Hubble slicing, we have 

S = 0(e 2 ). (A4) 

In general, we have 

i/> = O(l) . (A5) 

Also, from Eq. (|2.12|) . we have 

~ . a d t ip d t S d k tp k wp Q d k S k k k 

o t w + 6 — —it! H ; — w + 6 — —wv -I v + wo k v + (a k w)v — , (A6) 

if) T(l + 6) yj Tp 

where w = au°. Then Eq. (|Aip gives w = 1 + 0(e 2 ), and Eq. (|A6p gives 

d t i> = 0(e 2 ) . (A7) 

From Eq. (|2.10|) . we see 

Si, ~ 1 fS k k = 0(e 2 ) . (A8) 
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Therefore, from Eq. (|2 . 20[) we have 

dtAij = aKA tJ + <3(e 2 ). (A9) 
The uniform Hubble slicing condition (|2.23p . together with Eqs. (|A5[) and (|A7[) . gives 

X -a-l = 0(e 2 ). (A10) 

Then, Eq. (TA"9|) becomes 

dtAij = -3^ +0(e 2 ). (All) 

The lowest-order homogeneous solution of this equation is Aij oc a -3 , which could be 0(e). As discussed in the text, 
however, we assume 

dtJij = 0{e 2 ). (A12) 
This assumption, together with Eq. (|2.18p . implies that the 0(e) part of Aij is absent. We therefore have 

Ay = 0(e 2 ) . (A13) 

Finally, from Eq. (|2.16p . we have 

Uj =«°(« i +A) = 0(e 3 ), (A14) 
and inserting this into the expression for vP given in Eq. (|2.7[) gives 

w = au° = 1 + 0(e 6 ). (A15) 
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